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Starting from the Boltzmann equation in the relax;ation time approximation and employing a 
Chapman-Enskog like expansion for the distribution function close to eqnilibrium, we derive second- 
order evolution equations for the shear stress tensor and the dissipative charge current for a system 
of massless quarks and gluons. The transport coefficients are obtained exactly using quantum 
statistics for the phase space distribution functions at non-zero chemical potential. We show that, 
within the relaxation time approximation, the second-order evolution equations for the shear stress 
tensor and the dissipative charge current can be decoupled. We find that, for large values of 
the ratio of chemical potential to temperature, the charge conductivity is small compared to the 
coefficient of shear viscosity. Moreover, we show that in the relaxation-time approximation, the 
limiting behaviour of the ratio of heat conductivity to shear viscosity is qualitatively similar to that 
obtained for a strongly coupled conformal plasma. 


I. INTRODUCTION 

High-energy heavy ion collisions at the BNL Relativis¬ 
tic Heavy Ion Collider (RHIC) [1, 2] and the CERN Large 
Hadron Collider (LHC) [3-5] create strongly interact¬ 
ing matter under extreme conditions of high tempera¬ 
ture and density as it is believed to have existed in the 
very early universe [6, 7]. At such conditions, quarks 
and gluons are deconfined to form a new state of mat¬ 
ter, the quark-gluon plasma (QGP). The QGP behaves 
as a strongly coupled plasma having the smallest shear 
viscosity-to-entropy density ratio, ry/s [8-13]. Relativis¬ 
tic hydrodynamics has been applied quite successfully to 
describe the space-time evolution of the QGP formed in 
high-energy heavy ion collisions and estimate its trans¬ 
port coefficients [14]. 

In applications of hydrodynamics it is rather straight¬ 
forward to employ the ideal (Euler) equations. The in¬ 
clusion of dissipative effects in the evolution of the QGP 
started only a few years ago. However, most of the stud¬ 
ies have focused on exploring the effects of the shear 
viscosity on the QGP evolution and extracting its mag¬ 
nitude from experimental measurements. Nevertheless, 
there are other sources of dissipation such as bulk vis¬ 
cous pressure and dissipative charge current that may 
have a significant effect on the hydrodynamic evolution 
of the QGP. While the effects of bulk viscous pressure 
has been studied in some details [15-20], the dissipative 
charge current has been largely ignored. This may be 
attributed to the fact that at very high energies, baryon 
number and its corresponding chemical potential are neg¬ 
ligible. However, at lower collision energies such as those 
probed in the RHIG low-energy scan or at the upcoming 
experiments at the Facility for Antiproton and Ion Re¬ 
search (FAIR), baryon number can no longer be ignored 
and therefore charge diffusion may play an important 
role. 

The earliest theoretical formulation of relativistic dissi¬ 
pative hydrodynamics are due to Eckart [21] and Landau- 


Lifshitz [22]. However these formulations, collectively 
called relativistic Navier-Stokes theory, involve only first- 
order gradients and suffer from acausality and numerical 
instability due to the parabolic nature of the equations. 
Second order or extended theories by Grad [23], Muller 
[24] and Israel and Stewart (IS) [25] were introduced to 
restore causality. Therefore it is imperative that second 
order dissipative hydrodynamic equations should be em¬ 
ployed in order to correctly describe the evolution of the 
QGP. However, the IS formulation of a causal theory of 
relativistic hydrodynamics from kinetic theory, contains 
several inconsistencies and approximations, the resolu¬ 
tion of which is currently an active research area [26-38] . 

In order to formulate a causal theory of relativistic 
dissipative hydrodynamics from kinetic theory, it is de¬ 
sirable to first specify the form of the non-equilibrium 
phase-space distribution function. For a system close 
to local thermodynamic equilibrium, the non-equilibrium 
corrections to the distribution function can be obtained 
using either (i) Grad’s moment method [23] or (ii) the 
Ghapman-Enskog method [39]. Although both meth¬ 
ods involve expanding the distribution function around 
its equilibrium value, it has been demonstrated that the 
Ghapman-Enskog method in the relaxation time approx¬ 
imation results in a better agreement with microscopic 
Boltzmann simulations [32, 33] as well as with exact so¬ 
lutions of the Boltzmann equation in the relaxation-time 
approximation [32-36]. 

In the absence of conserved charges, the Ghapman- 
Enskog method has been used to compute the second- 
order transport coefficients for vanishing [32-34] as well 
as finite particle masses [35, 36]. On the other hand, in 
the presence of conserved charges but for classical par¬ 
ticles with vanishing masses, the second-order transport 
coefficients corresponding to charge diffusion (or alterna¬ 
tively heat conduction) have been obtained by employing 
the moment method [40, 41]. However, they still remain 
to be determined for quantum statistics. Here, we em¬ 
ploy the Ghapman-Enskog method to achieve this. 
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In this Letter, we present the derivation of second- 
order evolution equations for shear stress tensor and dis¬ 
sipative charge current for a system consisting of mass¬ 
less quarks and gluons. In order to obtain the form of 
the non-equilibrium distribution function, we employ a 
Chapman-Enskog like expansion to iteratively solve the 
Boltzmann equation in the relaxation time approxima¬ 
tion [32]. Using this expansion, we derive the hrst-order 
constitutive relations and subsequently the second-order 
evolution equations for the dissipative quantities. The 
transport coefficients are obtained exactly using quan¬ 
tum statistics for the quark and gluon phase-space dis¬ 
tribution functions with a non-vanishing quark chemical 
potential. Moreover, we show that, up to second-order in 
the gradient expansion, the evolution equations for the 
shear stress tensor and the dissipative charge current can 
be decoupled. We also find that, for large values of the 
ratio of chemical potential to temperature, the charge 
conductivity is small compared to the coefficient of shear 
viscosity. Finally we demonstrate that the limiting be¬ 
haviour of the heat conductivity to shear viscosity ratio, 
obtained here in the relaxation-time approximation, is 
qualitatively identical to that of a conformal fluid in the 
strong coupling limit. 


II. RELATIVISTIC HYDRODYNAMICS 


In the case of massless partons, i.e., massless quarks 
and gluons, the conserved energy-momentum tensor and 
the net-quark current can be expressed in terms of the 
single particle phase-space distribution function as [42] 

T^’' =j^P P^P'' idqifq + fq) + 9gfg\ 

= -h TT^^, (1) 

=Jdp [gqifq - fq)] = -f 71^, (2) 

where dp = (ip/[(27r)^jpj], is the particle four mo¬ 
menta, and Pq and Pg are the quark and gluon degeneracy 
factor, respectively. Here fq, fq, and fg are the phase- 
space distribution functions for quarks, anti-quarks, and 
gluons. In the tensor decompositions, e, P, and n are 
the energy density, pressure, and the net quark number 
density. The projection operator = p^^'^ — u^u'^ is 
orthogonal to the hydrodynamic four-velocity defined 
in the Landau frame: . We work with the 

Minkowskian metric tensor = diag(-|-, —, —, —). 

The dissipative quantities in Eqs. (1) and (2) are the 
shear stress tensor and the particle diffusion cur¬ 
rent . With the definition of the energy-momentum 
tensor in Eq. (1), the bulk viscous pressure vanishes 
in the massless case. The energy-momentum conserva¬ 
tion, = 0, and particle four-current conservation, 

dgN^ = 0, yields the fundamental evolution equations 


for e, and n, as 


e + (e + P)d ~ — 0, (3) 

(e -f P)w“ - V“P -k A“5^7r^‘^ = 0, (4) 

h + n9 + dfj,n^ = 0. (5) 


Here we use the standard notation A = u^d^A for co¬ 
moving derivatives, 9 = for the expansion scalar, 

-k - \9/\^^'' for the velocity stress 

tensor, and V“ = for space-like derivatives. 

In the following, we briefly outline the thermody¬ 
namic properties of a QGP in equilibrium. In this case, 
the phase-space distribution functions for quarks, anti¬ 
quarks and gluons are given by 


f{0) _ 1 

exp(,d u • p — a) -k 1 ’ 

(6) 

t(0) _ ^ 

exp{/3 u- p + a) + P 

(7) 

f(0) 1 

® exp(/3 u • p) — 1 ’ 

(8) 


respectively, where u ■ p = Ugp^, /3 = 1/T is the in¬ 
verse temperature and a = p,/T is the ratio of the quark 
chemical potential to temperature. We consider vanish¬ 
ing chemical potential for gluons because they are uncon¬ 
strained by the conservation laws. 

The temperature, T, and chemical potential, p, of the 
system is determined by the matching condition e = cq 
and n = uq, where eg and Uq is the energy density and 
the net quark number density in equilibrium. The energy 
density, pressure and the net quark number density for 
a system of massless quarks and gluons in equilibrium is 
given by 


eg = UgU,, Jdppf^p’' pq (4°) -k 4 °^) + 

(4<7g + 759 ) 77^^4 , Pqrr.2..2 , 9q .,4 
120 ^ + 4 ^ ^ 8^2^ 

gy f dpp^p'' Pq ( 4 °^ + 4 °’) + 9gfg°^ 


(9) 


Pn = -^A 


(45g-k 7g,)7r2 4 5g^2„2 


360 


12 ^ 




247r2 


no = UgJdpp^ Pq (^ 4 °^ - 4 °') 


6 ^ ^“^6772^ ■ 

The equilibrium entropy density then becomes 


( 10 ) 


( 11 ) 


^ eo + Po-pno ^ (4gg + 7gg)7r2 ^3 ^ 

T 90 6 ^ ^ 

The above expressions for eg, Pg, ng, and sg can also be 
obtained directly from the partition function of an ideal 
QGP [42], 


In Z = 


V 

T 


(4gg -k Tpq)'^'^ ^rp2^2 _|_ 9q ^4 


360 


12^ 


247r2' 


( 13 ) 
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where V is the volume of the system. Indeed, using the 
thermodynamic relations 


eo = 


no = 


dhiZ 

Y dT 

TdlnZ 

V dfi ’ 


T 

+ /rno, Po = —\nZ, 


So = 


V 

1 d{TlnZ) 

V ar ^ 


(14) 


one recovers Eqs. (9)-{12). The matching conditions 
e = eg and n = ng allows us to define thermodynamic 
quantities like temperature and chemical potential of a 
dissipative system. The pressure P can then be obtained 
from the equation of state of the system. 

Even if the equation of state, relating e, P and n is pro¬ 
vided, Eqs. (3)-(5) are not closed unless the dissipative 
quantities and are specified. However, before we 
derive the evolution equations for the dissipative quanti¬ 
ties, we need to obtain expressions for the derivatives of 
/3 and a. Using Eqs. (3)-(5) and Eqs. (9)-(ll), we get 


p = ^e + oY), d = oY), 

e + P e-\-P^ 


(15) 

(16) 


where represents terms which are of second or 

higher order in derivatives. Since dissipative forces are 
caused by thermodynamic gradients present in a non¬ 
ideal system, and are at least linear in the gra¬ 
dient expansion. Note that while Eq. (15) is terminated 
at first-order (sufficient for the present work), Eq. (16) is 
exact. 

The QGP is a strongly coupled system and is conjec¬ 
tured to be formed close to local thermodynamic equilib¬ 
rium. Therefore, the phase-space distribution function 
can be split into equilibrium and non-equilibrium parts, 
/ = where \Sf\/f^°^ <C 1. Hence, from Eqs. (1) 

and (2), the shear stress tensor and the particle dif¬ 
fusion current can be expressed in terms of Sf as 


Jdpp°‘p^ [9q{Sfq -k Sfq) + PgSfg] , (17) 

n'" = ^^Jdpp°‘[gqi6fg-Sfg)], (18) 


where = |(A^A^-kA^A^)-iA'^^A„^ isatraceless 
symmetric projection operator orthogonal to and A^i,. 
In the following, we obtain Sf up to first order by using 
the iterative solution of the Boltzmann equation in the 
relaxation-time approximation and then derive second- 
order evolution equations for the dissipative quantities. 


III. DISSIPATIVE EVOLUTION EQUATIONS 

The determination of the form of the non-equilibrium 
phase-space distribution function is a central problem in 
statistical physics. This can be achieved by solving a ki¬ 
netic equation like the Boltzmann equation. The Boltz¬ 
mann equation governs the evolution of the distribution 


function which provides a complete description of the mi¬ 
croscopic dynamics of a system in the dilute limit. The 
relativistic Boltzmann equation with the collision term 
written in the relaxation-time approximation is [43], 

P^dgf ={^f - , (19) 

where is the relaxation time. Note, that for different 
species of particles, with inter- and intra-species inter¬ 
actions, the relaxation-times are usually distinct. Thus, 
in general, one should consider the QGP as a true mul¬ 
ticomponent system. In the following, we consider the 
special case with a common relaxation time for all parti¬ 
cle species in the QGP. The general case with a different 
relaxation time is left for future work. 

We employ iterative solution of the Boltzmann equa¬ 
tion (19) to derive the dissipative equations [31-33]. The 
first-order expressions for shear stress tensor and dissi¬ 
pative charge current are obtained as [32], 

= =/3„T«V^a. (20) 


Here and /3„ are the first-order transport coefficients 
obtained after performing the momentum integrations in 
Eqs. (17) and (18). For a system of massless partons, as 
considered here, we find 


Y = 


e + P 


Jtg nY 

3 e + P’ 


( 21 ) 


where. 


+0 =3qj dp(u-p) 

6 27r^ a (tt^ -|- a^) 


( 22 ) 


Here = 1 — where r = 1 for Fermions (quarks 
and anti-quarks) and r = — 1 for Bosons (gluons). 

Using Eq. (20), we also obtain the first-order dissipa¬ 
tive corrections to the distribution function. 


r(0)?(0) 

Jq Jq 

Ao)+o) 
J q J q 

Ao)f(o) 
Jg Jg 


2{u-p)j5, 


-p^^p'^TTgu + 


2{u-p)l3, 

P 

2{u-p)p, 


-p^^p'^TTgu + 


-P^p'^ngy. 


n 


e+P 

u-pj 

n 

+ M 

e+P 

U-pJ 




(23) 


P^rif,, 


(24) 

(25) 


Note, that the tensorial form of dissipative corrections to 
the distribution function, as given in the above equations, 
is analogous to that of Grad’s 14-moment approximation 
[23]. However, the coefficients of these terms are differ¬ 
ent which has interesting implications in the context of 
relativistic heavy-ion collisions [38]. 
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To obtain the second-order evolution equations, we fol¬ 
low the procedure discussed in Ref. [26] . We consider the 
comoving derivative of Eqs. (17) and (18), and rewrite 
Eq. (19) in favour of 6f. Using Eqs. (23)-(25) and per¬ 
forming the momentum integrations, we finally obtain 
the second-order evolution equation for and n^. 




- Y 1 ^ ’ 


(26) 










3/3n 

e + P 




(27) 


Here — Vu^)/2 is the anti-symmetric 

vorticity tensor and we have ignored terms higher than 
quadratic order in the gradients [32]. Note that in 
the relaxation-time approximation, the Boltzmann relax¬ 
ation time Tn is the time scale for evolution of both 
and n'', i.e., = Tn = t^j. By comparing the hrst-order 

equations, Eq. (20), with the relativistic Navier-Stokes 
equations for the dissipative quantities [22], 


(28) 


the dissipative relaxation times can be related to the first- 
order transport coefficients = rj/PT^ and t„ = Kn/Pn- 

It is interesting to note that terms like and 

do not appear in Eq. (27) because their coefficients 
vanish. Note also that, with /3^ = (e -|- P)/5 calculated 
using the appropriate equation of state, Eq. (26) is valid 
even in a hadronic phase dominated by massless pions. 

The last term in Eq. (27) couples the evolution of 
dissipative charge current with the shear stress tensor. 
This type of coupling leads to disagreement with trans¬ 
port results as shown in Ref. [40]. We observe, how¬ 
ever, that using Eq. (20), the last term in Eq. (27) is, 

up to second-order in the gradient expansion, equivalent 
to —{6/5)n^a‘'^. Thus, the evolution equation, Eq. (27), 
for the dissipative charge current becomes 

/ \ 71 ^ Q 

n(f^) + _ = (29) 

Tn 5 

This is the main result of the present work. The compact 
form of the above equation makes it straightforward for 
direct implementation in a viscous hydrodynamic code. 
Note that while it is possible to formally rewrite any 
second-order equation only in terms of gradients using 
the first-order expressions, Eq. (20), it does not usually 
imply decoupling of the dissipative evolution equations 
[44]. However, we have ensured that the second-order 
terms in Eq. (29) is product of a dissipative quantity and 
a gradient, as inherent in Eq. (26). 

It is important to note that the gradient expansion con¬ 
verges only for small deviations from equilibrium. This 
implies that the second-order scheme is, strictly speak¬ 
ing, justihed only if the deviations from the hrst-order 



FIG. 1: (Color online) The ratio of charge conductivity to 
shear viscosity scaled by the temperature, KnT/rj, for two fla¬ 
vor (solid line) and three flavor (dashed line) massless quarks, 
plotted against n/T. 


relations, Eq. (20), are of second order, or smaller. This 
is true, in general, to ensure the convergence of gradi¬ 
ent expansion in the formulation of dissipative hydro¬ 
dynamics. Hence, for consistency, the initial conditions 
for the second-order evolution equations, Eqs. (26) and 
(29), should be chosen such that the constitutive rela¬ 
tions, Eq. (20), are satished. Nevertheless, it was shown 
in Refs. [33-36] that the solutions of the second-order dis¬ 
sipative hydrodynamic equations rapidly converge to the 
exact solution of the Boltzmann equation, irrespective of 
the choices of initial conditions for the dissipative quan¬ 
tities. Thus, the use of Eq. (20), to modify the second- 
order terms in the evolution equation for the dissipative 
charge current, is tenable even if the initial conditions 
are chosen such that these relations are violated initially. 


IV. TRANSPORT COEFFICIENTS 

The dimensionless ratio /rj is a measure of the rela¬ 
tive importance of the charge conductivity and the shear 
viscosity. This quantity, in the relaxation-time approx¬ 
imation, is given by UnT/rj = /3„T//3.n., which can be 
studied as a function /r/T. To quantify this ratio, one 
still need to specify the appropriate quark and gluon de¬ 
generacy factors, Qq and Qg, as 

gq = Ns xNcxNf = 6Nf, 

gg = NsX {n! - 1) = 16, (30) 

where Ng = 2 is the number of spin degrees of freedom, 
Nc = 3 is the number of colours, and Nf is the number 
of flavours. 

In Fig. 1, we show the ratio /g as a function of fj,/T 
for Nf = 2 and Nf = 3. We observe that, while for small 
/r/T, this ratio is almost constant, it drops rapidly for 





5 



FIG. 2: (Color online) The ratio of thermal conductivity to 
shear viscosity, Kqfrj, scaled by the factor for two fla¬ 

vor (solid line) and three flavor (dashed line) massless quarks, 
plotted against fJ./T. 


statistics and it does not appear for a classical Boltzmann 
gas. In the limit of small /i/T, the constant C in Eq. (33) 
is C = {Agg + 7gq)l9gg. Thus, C = 37/27 and 95/81 
for two and three flavor QGP, respectively; see Fig. 2. 
On the other hand, for large g/T we find C = 5/3, in¬ 
dependent of the number of flavors, as shown in Fig. 2. 
These values of C are comparable to C = 8/9 (here a fac¬ 
tor of 1/9 indicates that the baryon chemical potential 
is three times the quark chemical potential) obtained in 
the calculations for strongly coupled conformal plasmas 
with finite chemical potential [45]. However, it should 
be noted that for a strongly coupled conformal plasma, 
the coefficient C depends on the number of space-time 
dimensions and is shown to be equal to 32/9, 8/9 and 
2/9 for four, five and seven dimensions, respectively [47]. 
Nevertheless, it is intriguing that up to a constant of 
proportionality, the limiting behaviour of the ratio Kq/g^ 
obtained in the relaxation-time approximation, is exactly 
the same as that derived in the case of a strongly coupled 
conformal fluid. 


larger g,/T, indicating that the conductivity of the QGP 
is small relative to the shear viscosity at low temperature 
and high density. Although the qualitative behaviour 
remains the same for Nf — 2 and Nf = 3, the drop 
in KnT/g is more pronounced for Nf = 3. Moreover, 
we note that at fi/T > 1 the ratio KnT/g is almost Nf 
independent. 

A further interesting quantity is the ratio of thermal 
conductivity to shear viscosity. The heat flow is related 
to the dissipative charge current via the relation g^ = 
— (e -I- P)nf^/n, and is given as [45] 

nT^ ^ e + P\^ 

= -Kq ^ Kq = Kn {-^ j > (^l) 

where Kq is the coefficient of thermal conductivity. Using 
the first-order relations in the relaxation-time approxi¬ 
mation, we find 


Kq _ (3n / £ + P\ 
g /3^ V nT / 


(32) 


In Fig. 2, we plot Kq/g scaled by the factor g?/' k'^T ver¬ 
sus g/T, for a two and three flavor QGP. We observe a 
constant behaviour in the limit of small as well as large 
g/T. Moreover, for large g/T, we see that {Kq/g)g‘^/n'^T 
is independent of number of flavors. These limiting be¬ 
haviours have interesting consequences. 

In the limit of both small and large g/T, Eq. (32) 
reduces to 


Kq 

1 



(33) 


which is similar to the Wiedemann-Franz law [45, 46]. 
The factor in the above equation is due to quantum 


V. CONCLUSIONS AND OUTLOOK 

In this paper we employed the iterative Ghapman- 
Enskog method to derive the second-order dissipative hy- 
drodynamical equations for a system of massless quarks 
and gluons. The bulk viscous pressure vanishes for such a 
system and therefore the dissipation is solely due to the 
shear stress tensor and the dissipative charge current. 
For the equilibrium distribution function, we considered 
quantum statistics with non vanishing quark chemical 
potential. We obtained novel, exact relations for the 
second-order transport coefficients corresponding to the 
dissipative charge current evolution. 

Moreover, we demonstrated that the evolution equa¬ 
tions for shear stress tensor and dissipative charge cur¬ 
rent can be decoupled. We also found that, for large 
values of the ratio of chemical potential to temperature, 
the charge conductivity is small relative to the shear vis¬ 
cosity. Finally, we showed that in the relaxation-time ap¬ 
proximation, the limiting behaviour of the ratio of heat 
conductivity to shear viscosity is qualitatively similar to 
that of a conformal fluid in the strong coupling regime. 

At this juncture, we would like to stress that the itera¬ 
tive Ghapman-Enskog approach employed here to obtain 
the dissipative evolution equations from the Boltzmann 
equation in the relaxation-time approximation is compat¬ 
ible with the gradient expansion inherent in the formula¬ 
tion of dissipative hydrodynamics, as opposed to the mo¬ 
ment method [32]. Looking forward, it would be interest¬ 
ing to determine the effect of the dissipative charge cur¬ 
rent in high-energy heavy-ion collisions, by implement¬ 
ing the dissipative equations derived here, in realistic hy¬ 
drodynamic simulations. A further challenging problem 
would be to extend the current second-order formulation 
to third order. We leave these questions for future stud¬ 
ies. 
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